Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 4 February 2008 (MN MfeK style file vl.4) 



Gravitational instability of discs with dissipative corona 
around supermassive black holes 

Fazeleh Khajenabi 1 * and Mohsen Shadmehri 2 ' 3 t 

1 School of Mathematical Sciences, University College Dublin, Belfield, Dublin 4, Ireland 

2 School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland 

3 Department of Physics, School of Science, Ferdowsi University, Mashhad, Iran 

o 

O ' 4 February 2008 

(N 

ABSTRACT 

We study the dynamical structure of a self-gravitating disc with corona around a 
super-massive black hole. Assuming that the magneto-rotational-instability (MRI) 
' responsible for generating the turbulent stresses inside the disc is also the source for a 

magnetically dominated corona, a fraction of the power released when the disc matter 
accretes is transported to and dissipated in the corona. This has major effect on the 
structure of the disc and its gravitational (in)stability according to our analytical and 
self-consistent solutions. We determine the radius where the disc crosses the inner 
radius of gravitational instability and forms the first stars. Not only the location of 
£vq , this radius which may extend to very large distances from the central black hole, but 

£T) ■ also the mass of the first stars highly depends on the input parameters, notably the 

f^) viscous coefficient, mass of the central object and the accretion rate. For accretion 

t~-» , discs around quasi-stellar objects (QSOs) and the Galactic center, we determine the 

■ self-gravitating radius and the mass of the first clumps. Comparing the cases with a 

"*^h \ corona and without a corona for typical discs around QSOs or the Galactic center, 

O { when the viscosity coefficient is around 0.3, we show that the self-gravitating radius 

I 1 decreases by a factor of approximately 2, but the mass of the fragments increases with 

more or less the same factor. Existence of a corona implies a more gravitationally 
unstable disc according to our results. The effect of a corona on the instability of the 
disc is more effective when the viscosity coefficient increases. 
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1 INTRODUCTION 2007). Quasar spectra show that near the black hole there is 

high metallicity (Hamann et al. 2002), presumably because 

Among various physical agents which play significant roles r , „ ,. T , ,. , , 

° . or massive star formation, interestingly, young and massive 

in accreting systems, self-gravity has a vital role in accre- , , ., , . „ , r „ . * . „ , , . 

, ° stars are observed within 0.5 pc of Sgr A in our Galactic 

tion discs around supermassive black holes. According to , T , ,-. , . , „„„„ ,-, i 

or, center (e.g., Lu et al. 2005; Paumard et al. 2006; Bender et 

the standard theory of accretion discs Shakura & bunyaev , nr ,„ r \ n , , .. c ^ , .. 

J v J al 2005). Based on the observation of the Galactic center, 

1973), the outer parts of steady, geometrically thin, opti- T . „ „ , , , , nrvrvo \ , , . , , c 

„ , , , , , , Levin & Befoborodov (2003) proposed that a recent burst of 

calfy thick discs are prone to seff-gravity and they might r , . , , . , , . , 

' . „, o t. , s t ar formation has occurred m a dense gaseous disc around 

be expected to fragment into stars (e.g., Shlosman & Begel- ggr A * ^ A&t& q{ p aumard et al (2006) strongly 

man 1987; Goodman 2003). While many authors proposed r ., r ,. „ , ,. . , , 

' . , . , , . „ „„ils . favour m situ star formation from dense disc in the centraf 



regions of the Galactic center. Star formation activities are 
seen in the nuclei of Seyfert galaxies (Terlevich 1996). 



possible solutions to this problem (e.g., Goodman 2003), it 
seems such a gravitationally unstable disc is a good expla- 
nation for the existence of the first starts in galactic centers 
or the formation of supermassive stars in quasar accretion Thege observations have mot ivated many authors to 

discs (Goodman & Tan 2004; Tan & Blackman 2005; Levin gtudy gtar formation under extrem e conditions in discs 

around massive black holes at the center of galaxies. The 
suggestion of star formation in the self-gravitating parts of 
* E-mail: fazeleh.khajenabi@ucd.ie (FK); an accretion disc was first made by Kolykhalov & Sunayev 

t E-mail: mohsen.shadmehri@dcu.ie (MS) (1980). Shlosman & Begelman (1989) discussed in detail the 
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conditions for star formation in the discs around black holes. 
Recently, star formation near to the Galactic Center has also 
been studied by Levin (2006) and Nayakshin (2006). Levin 
(2006) considers a disc that fragments into many stars of 
up to a few hundred solar masses and Nayakshin (2006) 
shows that the initial mass function for disc-born stars at the 
distances between 0.03 and 0.3 parsec from the supermas- 
sive black hole should be top-heavy. Collin & Zahn (1999a, 
1999b) have also studied formation of stars in AGN discs. 
Goodman & Tan (2004) proposed that supermassive stars 
may form in quasar discs. According to their model, while 
the disc may fragment and each fragment may start with 
only a hundred solar mass, it seems likely they grow to a sig- 
nificant fraction of the disc mass, perhaps around 10° solar 
mass. Just recently Nayakshin et al. (2007) studied star for- 
mation in the central parsec of our Galaxy using numerical 
simulations. However, they found that the masses of the first 
clumps are much smaller than 10 5 solar mass which has been 
predicted by simple theoretical analysis (e.g., Goodman & 
Tan 2004). 

The observation of the UV bump of AGNs (e.g., Malkan 
& Sargent 1982) argues in favour of geometrically thin and 
optically thick discs, possibly embedded in a hot X-ray emit- 
ting corona. In fact, coexisting soft and hard X-ray compo- 
nents observed in ordinary AGNs can not be well described 
by simple one-zone hot-flow models and require composite 
disc-corona structure. Considering some common features of 
the corona of a disc with a solar corona, interesting theoret- 
ical investigations have also been presented in the frame of 
a magnetized disc corona (Haardt &: Maraschi 1991; Field 
& Rogers 1993; Nakamura & Osaki 1993; Merloni & Fabian 
2002; Misra & Taam 2002; Liu, Mineshige & Shibata 2002; 
Merloni 2003) . It has been suggested that the corona consists 
of localized active regions. It is likely that these are produced 
by magnetic fields in the disc amplified through differential 
rotation. When the disc's magnetic field builds up signif- 
icantly, buoyancy forces the field out and above the disc, 
giving rise to active regions of high magnetic field. Misra & 
Taam (2002) found that a significant part of the accretion 
flow (or dissipation rate) can take place in the corona if the 
scale height of the magnetic field is larger than that of the 
disc. According to their study, optically thick discs with dis- 
sipative coronas can provide an attractive explanations for 
the origin of the soft spectral component observed in black 
hole X-ray binary systems. Rozariska et al. (1999) studied 
the vertical structure of a radiation-pressure-dominated disc 
with a hot corona. They showed that the presence of the 
corona modifies considerably the density and the opacity of 
the disc surface layers and concluded that the corona is not 
only essential from the point of view of the X-ray spectra 
formation but helps to remove the problems with the disc 
models for AGN. The disc-coronal model for AGN has also 
been discussed by Field & Rogers (1993). 

Merloni & Fabian (2002) showed that, if angular mo- 
mentum transport in the disc is due to magnetic turbulent 
stresses, the magnetic energy density and effective viscous 
stresses inside the disc are proportional to the geometric 
mean of the gas pressure and the total pressure (gas plus 
radiation) . They discussed why energetically dominant coro- 
nae are ideal sites for launching powerful jets/outflows, both 



MHD- and thermally driven. Merloni (2003) presented thin 
disc solutions accompanyied by powerful, magnetically dom- 
inated coronae and outflows, as models for black holes ac- 
creting at super-Eddington rates. Based on the assumption 
that the magnetorotational instability (MRI) responsible for 
generating the turbulent stresses inside the discs (Balbus & 
Hawley 1991) is also the source for a magnetically dominated 
corona, Merloni & Nayakshin (2006) studied the limit-cycle 
instability in magnetized accretion discs. In these models, 
a fraction of the gravitational energy release is assumed to 
dissipate in a hot diffuse corona, above the main body of the 
disc, away from the majority of the accreted mass. 

Thus, one may conclude that corona is an important 
part of the accretion processes around central black holes 
of AGNs. Possible effects of a corona on its underlying disc 
have been studied using different models. For example, a 
dissipative corona may be able to thermally stabilize the 
underlying disc (Misra & Taam 2002). However, the possi- 
ble effects of the corona on the gravitational stability of the 
disc has not been studied to our knowledge. In this study, we 
would like to address this question by a very simple model: 
What is the effect of the corona on the gravitational sta- 
bility of a self-gravitating disc? We know that the corona 
may thermally stabilises the disc, but can we expect gravita- 
tional stabilization as a result of the presence of the corona? 
Considering the growing interest in star formation near our 
Galactic center and discs of AGNs, any possible effects of the 
corona will have interesting implications in these studies. In 
our approach, the disc structure equations are reduced to 
simple algebraic relations in terms of average quantities of 
the disc and corona. 



2 GENERAL FORMULATION 

We assume that all angular momentum transport takes place 
in the disc and the mass accretion rate M is constant with 
radius and time. However, the microphysics mechanism of 
the angular momentum transfer remains unknown. Shakura 
& Sunyaev (1973) replaced all the missing physics by a pa- 
rameter a. This approach has been widely used for study- 
ing the dynamics and structure of the accretion flows. A 
promising mechanism for driving the turbulence responsible 
for angular momentum and energy transport is the action 
of the magneto-rotational instability (MRI) that is expected 
to take place in such discs (Balbus & Hawley 1991). But for 
a disc-corona system it is typically thought that the viscous 
stress, assumed to be magnetic in nature, transports angular 
momentum and initially randomise the gravitational binding 
energy near the midplane. The magnetized fluid elements, 
which are buoyant with respect to their surroundings, dissi- 
pate above the disc. Merloni & Fabian (2002) and Merloni 
(2003) presented very detailed discussions about dissipation 
in the corona and its relation with angular momentum trans- 
port in the disc itself. Our disc-corona model is developed 
along the line proposed by Merloni & Fabian (2002) and 
Merloni (2003). 

We consider a more general prescription for the viscous 
stresses T T( p (Taam & Lin 1984; Watarai & Mineshige 2003; 
Merloni & Nayakshin 2006): 
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1 — till u/2 



(1) 



where ao and < /i < 2 are constants. Also, p is the sum of 
the gas and radiation pressures. Phenomenological models 
generally assume that at each radius, a fraction / of accre- 
tion energy is released in the reconnected magnetic corona. 
Assuming that in MRI-turbulence discs such a fraction / of 
the binding energy is transported from large to small depths 
by Poynting flux, Merloni & Nayakshin (2006) estimated the 
fraction / as 



(2) 



where /3 is the ratio of gas pressure to the total pressure. So, 
in this model, this fraction / is not a free parameter. 

The rotation curve is dominated by a Newtonian point 
mass M, as relativistic effects are only important at small 
radii. Although self-gravity may be important for local sta- 
bility at the outer radii of the disc, the rotation angular 
velocity is not much affected because the disc is thin and its 
mass is smaller than the central point mass (e.g., Goodman 
2003). Thus, the rotational angular velocity of the disc is 
Keplerian, i.e. S!k = \J GM/R 3 . Having equation (1) as a 
prescription for the viscous stresses and equation (2) as an 
expression for the fraction of power dissipated in the corona, 
we can write the basic equations of the disc. Vertical hydro- 
static equilibrium of the disc implies 



£ 



(3) 



and the azimuthal component of the equation of motion 



3n K MJ(R), 



(4) 



where J(R) = 1 — \J R ln /R- Here, R ln denotes the inner 
boundary of disc. Since we are interested in the regions of 
the disc with radii much larger than Ri n , we can assume 
J « 1. However, we will keep this factor in order to present 
solutions of the disc structure as general as possible. 
The energy equations is given by 



ar c 4 ff = ffi|MJ(i?)(l-/), 



(5) 



where as we discussed / = y // 2aof3' J -/ 2 . Also, we can assume 
that the vertical transport of heat is by radiative diffusion 
which implies the midplane and the surface temperatures 
are related by 



r = (^E) 1 / 4 Teff , 



(6) 



where k is the opacity coefficient. 

Now, equations (3), (4) and (5) are the main equations 
which enable us to find p and T as functions of R and /3 and 
the other input parameters. Thus, 



_ AoQk 1/2 16-7r 2 Q ( 2 )Cfc B 

3nao Safi^mnM^J 2 ' 



-1/3 



(1 - ^0^/2)1/2 )/l2 
(1-/9)1/3 1 



(7) 



167T 2 a^cks 



3k«o 3a/i m m H M 2 f24 J 2 



-2/3 



(1 - Vjogf; .(8-5 M )/12 
(1-/3)2/3 P 

There is an algebraic equation for f3 as follows 

167r 2 ao c ^B 



(8) 



^ m m H 3fta 



-5/3 



(1- y/2a (}»/*) 3 ' 2 p (8+ll)/12 



= 1, 



(9) 



(1-/3)5/3 

where /i m is the mean particle mass in units of the hydro- 
gen atom mass, m H . The other constants have their usual 
meanings. 

In order to study behavior of our solutions, it is more 
convenient to introduce dimensionless variables. For the cen- 
tral mass M, we introduce M 8 = M/(1O 8 M ) and for 
the radial distance R, we have — R/(10 i Rs), where 
Rs — 2GM/c 2 is the Schwarzschild radius. The mass ac- 
cretion rate can be written as 

IE Ib Lb 



M = 



(10) 



where Ie = L/Le is the dimensionless disc luminosity rela- 
tive to the Eddington limit and e = L/(Mc 2 ) is the radiative 
efficiency. Also, k c . b . w 0.04 m 2 kg _1 is the electron opac- 
ity. In our analysis, we will use the nondimensional factor 
Ze/c as a free parameter so that by changing this parame- 
ter we can consider appropriate values of the accretion rate. 
However, some authors introduce different forms for the ac- 
cretion rate. For example, Nayakshin & Cuadra (2005) who 
studied gravitational stability of the Galactic center, intro- 
duced M = (m/e)(L E /c 2 ) with e w 0.06 and rh = 0.03 to 
1. These values, which are appropriate for the Galactic cen- 
ter, correspond to Ze/c ~ 0.5 to 16.6 in our notation. Also, 
for a central mass with mass M = 10 s Mq, Goodman & Tan 
(2004) proposed Ze/e = 10. Thus, in our analysis, the chosen 
values of Ie/s = 1 and 10 are acceptable. 

Now, we can rewrite our solutions as (in SI) 

p = 2.76 x 10- 6 «- 1/2 ^ /2 M 8 - 1/2 (^) 2 J 2 r 3 - 15/4 



x/3 (8-M)/4 (1 _ py2 {1 _ ^2^/2)3/2 > 

p = 27.15a - 5/6 K 1/2 M- 5/6 (^) 4/a J 4/3 r- 13/4 



(11) 



(12) 



(13) 



x/3 (8-5 M )/12 (1 _ ^-2/3^ _ y / 2a0 P^) 1 ' 2 , 

f = 4.67 x 10- 3 ao - 1/6 A- 1/2 M 8 - 1/6 (V 1/3 J~ 1/3 

x/r (8+M)/12 (1 _ ^2/3^ _ v / 2qo/3p /2 ) -1/2 ] 

and the ratio (3 is obtained from nondimensional form of 
equation (9), i.e. 

0.16« ( 1 /^ 3/2 M 8 1/6 (^) 4 / 3 J 4 / 3 r 3 - 7/4 

x/3 (8+m)/12 (1 _ ^"5/3^ _ ^ 2ao/3 M/2)3/2 = ^ 



(14) 



where k = k/k c . s . and we assumed /i m « 0.6. We can also 
calculate the surface density as 



5„-2/3~ Aif l/3^5/3 7 5/3 r -ll/4 



J .27 x 10° a""' ^kM^^-jY'-l 
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x(3 



(4-A0/3 



(l-/3)- 4/3 (l- VW^). 



(15) 



Equations (11), (12) and (13) along with equation (14) 
describe the structure of a disc with a dissipative corona. 
However, the physical variables depend not only on the ra- 
dial distance but also on the ratio of the gas pressure to 
the total pressure which can be calculated at each radius 
from algebraic equation (14). For (j, = 1, our solutions re- 
duce to what has been obtained by Merloni (2003). In the 
next section we will analyze our solution, in particular the 
gravitational stability of the disc. 
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3 ANALYSIS 

Since we are interested in gravitational stability of the disc, 
we can approximate J(R) — 1 for R ^> R ln . Having the 
input parameters, we can solve equation (14) numerically 
at every radius to obtain the ratio of the gas pressure to 
the total pressure. Then, equation (2) gives the fraction 
/ of gravitational power associated with the angular mo- 
mentum transport that is transported vertically and dis- 
sipated in the corona. As the systems tends toward the 
gas pressure-dominated regime, the dissipated energy in the 
corona increases according to this equation. However, our 
solutions generally describe an inner radiation-dominated 
region with an outer gas pressure-dominated region accord- 
ing to equation (14). Moreover, the size of gas pressure- 
dominated regime increases, as the accretion rate decreases. 
These results are valid irrespective of the value of \i the ex- 
ponent of the magnetic viscosity. For fj, = 1, similar typical 
behaviors have already been obtained by Merloni & Fabian 
(2002). We want to extend this analysis by studying the 
gravitational instability of the disc. In doing so, we should 
calculate the Toomre parameter of the disc. 

Toomre (1964) showed that a rotating disc is subject to 
gravitational instabilities when the Q-parameter becomes 
smaller than a critical value, which is close to unity, 



Q 



ttGE' 



(16) 



where c s is the sound speed inside the accretion disk and 
Q = Qk is the angular velocity. So, the Toomre parameter 
of our model becomes 



Q = Ua 



1/2.-3/2 



x/r 



(I 



-M)/4 



ill. 



-3/2/ 



(-)-V 



-2 3/2 



(1-/3) 2 (1- V2ao^)" 3/2 



(17) 



This equation with algebraic equation (14) gives the Toomre 
parameter as a function of the radial distance. Generally, 
this parameter is much higher than unity in the inner parts 
of the disc which implies these regions are gravitationally 
stable and do not fragment. But the Toomre parameter de- 
creases with increasing radial distance so that Q reaches the 
critical value of unity at a self-gravitating radius which we 
denote by _R sg . Thus, all regions with R > _R sg are gravita- 
tional unstable and may fragment to clumps and cores. 

Different authors estimate the mass of fragments differ- 
ently. Since the disc is marginally unstable, the initial sizes 
and masses of gravitationally bound fragments can be de- 
termined by Toomre's dynamical instability (Toomre 1964) . 
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Figure 1. The self-gravitating radius R Bg (in Schwarzschild ra- 
dius ijg) and the mass of the first clumps Mf rag (in solar mass) at 
this radius vs. the exponent fi of the viscosity prescription for a 
central black hole with mass M = 10 8 Mq and k = 1. Each curve 
is marked by a pair (ao, 2eA) ■ 



The most unstable wavelength for the Q ~ 1 disc is of order 
of the disc vertical scale height H (Toomre 1964). Thus, the 
most unstable mode has radial wave number fc mu = (QH)~ X 
and so the mass of a fragment at R = i? sg becomes 



M frag » E(^) 2 = 47r 2 EH 2 . 



(18) 



3.1 Parameter study 

Now, we can have a parameter study of the solutions. Fig- 
ure 1 shows the self-gravitating radius i? sg (top) and the 
mass of the fragments (bottom) versus the exponent \x of 
the viscosity prescription for different input parameters. We 
fix the central mass (Ms = 1) and the opacity (R = 1), but 
vary the other input parameters. Each curve is labeled by 
the pair of the viscosity coefficient ao and the accretion rate 
Zs/e, as (ao, te/e). The top panel of Figure 1 shows that the 
disc of our system becomes self-gravitating at distances close 
to the central black hole. While for high accretion rate the 
self-gravitating radius R sg depends on the exponent fx, for 
lower accretion rate this radius is more or less independent 
of the exponent of the viscosity. For a fixed accretion rate, 
the self-gravitating radius R ss increases with ao. 

In fact, our viscosity prescription is different from the 
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Table 1. The self-gravitating radius R Bg (in Schwarzschild radius R%) and the mass of the first clumps 
Mf rag (in solar mass) for M = 10 s Mq, fi = 1 and k = 1. 



without corona 




We = 1 






= 10 






Q!0 


0.03 


0.3 


0.03 






0.3 




R S g/{10°Rs) 


0.5 


3.8 


0.75 






2.3 




M frag /M Q 


993 


17 


14350 






357 




with corona 




l E /e = 1 






= 10 






«0 


0.03 


0.3 


0.03 






0.3 




R sg /(10 3 i? s ) 


0.4 


1.7 


0.66 






1 




M frag /M 


1104 


29.7 


15530 






577 






Table 2. The same as Table 1, but for M = 3 X 10 6 M Q . 


without corona 




Ijs/e = 1 






ZbA = 


10 




no 


0.03 


0.3 




0.03 






0.3 


R s J(W i R s ) 


470 


6884 




102 






1484 


M frag /M 


0.004 


5 x 10- 


5 


0.12 






0.0016 


with corona 




l E /e = 1 






w« = 


10 




ao 


0.03 


0.3 




0.03 






0.3 


R sg /(10 3 R S ) 


407 


3275 




88 






704 


Affrag/M 


0.0045 


10" 4 




0.14 






0.003 



standard alpha model (Shakura & Sunyaev 1973). In alpha 
model, Rice et al 2005 argue that the maximum viscosity 
in marginally stable self-gravitating discs is around 0.06. 
Unfortunately, there is not enough information about ap- 
propriate values of coefficient ao for studying gravitational 
structure of discs. In analogy to previous studies of accretion 
discs based on the viscosity prescription of our model (e.g., 
Watarai & Mineshige 2003), we are using a range of 0.03 to 
0.3 for ao. 

Using equation (18), we can calculate the mass of the 
first fragments for the set of the above input parameters. For 
low accretion rate, the mass the fragments is also more or 
less independent of the exponent of fi. But low values of ao 
correspond to fragments with higher masses so that for pair 
of (0.03, 10) the masses will be between 10 4 and 10 5 solar 
mass depending on the viscosity exponent /x. However, as 
we mentioned, we think low values of ao are not acceptable 
in self-gravitating discs. For a fixed accretion rate, the mass 
of the first clumps decreases with increasing ao- But as the 
accretion rate decreases, the mass of clumps decreases as 
well so that while for (0.3, 10), we have Aff rag ~ 600Mq for 
(0.3, 1) the mass of the fragments decreases approximately 
to Mf rag ~ 3OAf0 . It means for the above set of input param- 
eters, as the accretion rate decreases, fragments with lower 
masses are forming at larger distances from the central black 
hole. According to our results not only the self-gravitating 




R/(10 3 RJ 



Figure 2. The integral fraction / of accretion power dissipated 
into the corona vs. radial location in the disc with fj, = 1.5 for 
M = 10 s Mq (solid lines) and M = 10 7 Mq (dashed lines) central 
black hole with £t m = 0.6 and k = 1. The curves are labeled by 
the ratio ZeA- A s this ratio increases, the fraction / dissipated 
energy in the corona decreases, irrespective of the other input 
parameters. 
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radius but also the mass of the first clumps show wide ranges 
of variations depending on the input parameters. In order 
to make easier comparison with other input parameters, Ta- 
ble 1 summarizes the self-gravitating radius and the mass of 
the first clumps for M = 10 8 M Q , [i = 1 and k = 1. We will 
discuss the fraction of dissipated energy in the corona. But 
before that we can study the solutions for other set of the 
input parameters. 

We found the self-gravitating radius and the mass of the 
fragments are not very sensitive to the exponent fi for lower 
mass of the central black hole. Thus, we summarize our re- 
sults for M = 3 x 10 6 M Q , /j, = 1 and k = 1 in Table 2 which is 
appropriate for modeling the Galactic center. Clearly, when 
the mass of the central black hole decreases, the disc is more 
gravitationally stable as R sg is much larger than in Figure 
1 and the mass of the fragments are much smaller. In this 
case, also the typical behaviors of the solutions are similar 
to the case with M — 1O 8 M0. We see that as the viscosity 
coefficient ao increases, not only does the self-gravitating ra- 
dius increase, but the mass of the fragments decrease. Also, 
as the accretion rate increases, the self-gravitating radius 
decreases and the mass of the fragments increases as well. 

Goodman & Tan (2003) also estimated the mass of first 
clumps in quasar discs. Clearly, the disc of their model is 
without corona. Assuming that the viscosity is proportional 
to gas pressure, they found R sg w 1700J? S to f? sg « 2700i? s 
with a few hundred solar mass for the fragments. Another 
point is that their model is radiation pressure dominated. 
In another study, Nayakshin (2006) studied star formation 
near to our Galactic center. The mass of the first stars of 
this model is a few solar mass which may increase because of 
the subsequent accretion. However, our model shows wider 
ranges for the self-gravitating radius and the mass of the 
fragments. Moreover, when we reach the self-gravitating ra- 
dius, the disc of our model is not necessarily radiation pres- 
sure dominated. In fact, we found the transition to the self- 
gravitating regime happens at radii where generally, the gas 
pressure, if it is not dominant, is at least comparable to the 
radiation pressure. 

We can calculate the fraction / of dissipated energy in 
the corona. Instead of integrating the fraction / over all the 
disc, for our purposes it is sufficient to plot this fraction as 
a function of the radial distance. Figure 2 shows the typi- 
cal behavior of / for our sets of the input parameters. Note 
that since we are not interested in the very inner regions of 
the disc, this plot is not accurate near these regions as we 
approximated J«l. This figure shows interesting behavior. 
First of all, as the viscosity coefficient increases, the fraction 
/ increases as well. We also showed the self-gravitating ra- 
dius increases and the mass of the fragments decreases with 
increasing ao- It means the disc is more gravitationally sta- 
ble as ao increases or correspondingly, the amount of the 
dissipated energy in the corona increases. In order to con- 
firm this result, we should check the behavior of / and the 
disc for the other input parameters. 

Consider Figure 2 for a fixed coefficient ao- We see 
that as the accretion rate decreases, the fraction / increases 
which means more energy dissipates in the corona. Also, as 
the mass of the central black hole decreases, the fraction 
/ increases for a fixed accretion rate. For these cases, we 



showed that the disc becomes self-gravitating at larger radii 
with lower masses for the fragments, as the mass of the cen- 
tral black hole decreases (compare with Figure 1 and Table 
!)• 

3.2 The effect of corona 

How can we understand the above results based on the ex- 
istence of the corona? For addressing this question, we can 
either compare our results directly with previous authors or 
with our calculations for the standard case (no corona). Al- 
though we briefly mentioned previous studies in this field, 
e.g. Goodman & Tan (2004) and Nayakshin (2006), we can 
not make direct comparisons with those studies because our 
basic assumptions for the corona and viscosity prescription 
are different. We think a straightforward way to understand 
the effect of the corona on the gravitational instability of the 
disc is direct comparison with our solutions, but without a 
corona. 

Obviously, a certain amount of the generated energy in 
the disc transports in the corona in our model. This leads to 
a cooler disc in comparison with a disc without corona. One 
may conclude that since the sound speed and consequently 
the Toomre parameter decrease, a disc with corona will be 
more unstable gravitationally. This result is correct as long 
as the surface density does not change or increase. In our 
model, the surface density is also lower than in the case of a 
disc without corona. So, this qualitative consideration does 
not help us to determine the typical behavior of the Toomre 
parameter as a result of existence of the corona. 

Now, we can simplify equation (17) by calculating rs 
from equation (14), 

Q = 9.14a?/ 14 £- 3/14 M 8 - 19/14 (^)- 6 / 7 

x/3 (9 M -40)/28 (1 _ p) 4/7 {1 _ ^/ 2ao f3^r 3/14 . (19) 

Also, the Toomre parameter Qb of the disc without a corona 
becomes 

Qb = 9.Ual /14 k- :i/14 M- 19/14 ( l -^)- 6/7 

x/3 (9m-40)/28 (1 _ /3) 4/7^ (2q) 

Having all the input parameters fixed at each (3 there is 
a relation between Toomre parameters of discs with (and 
without) corona: 

^(i-^r 3/i4 . (21) 

Qb 

Since the expression inside the above parenthesis is always 
lower than unity, the above equation implies Q > Qb for 
each (3. On the other hand, numerical solutions of equa- 
tion (14) for discs with (and without) corona show that 
7*3, b > 7*3 for a fixed /3. It means that the disc without corona 
reaches to a fixed (3 at larger radius in comparison with a 
disc with corona. When we plot the Toomre parameter ver- 
sus the distance from the central black hole, we have always 
Qb > Q at a fixed radius. Thus, the Toomre parameter 
Q of the disc with corona reaches unity at smaller radius 
in comparison with a disc without corona. In other words, 
the self-gravitating radius of a disc with corona is smaller 
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than disc without corona. We can confirm this general re- 
sult by direct numerical solution of the equations. Figure 
3 shows typical behavior of the Toomre parameter for the 
cases with and without a corona. Also, Tables 1 and 2 show 
the self-gravitating radius and the mass of the fragments of 
the discs with and without corona. In the last subsection we 
discussed general properties of discs with corona. Now, we 
can compare the gravitational instability of the discs with 
and without a corona. 

By comparing the Tables, we can simply conclude that 
as the mass of the central black hole increases, the discs 
become more self-gravitating irrespective of the existence 
of corona. In fact, while the self-gravitating radius i? sg de- 
creases with the mass of the central object, the estimated 
mass of the clumps increases. These Tables confirm our 
above analytical expectation that in discs with a corona the 
Toomre parameter is lower than the same disc but without 
a corona. For example, Table 1 compares self-gravitating ra- 
dius and the mass of the clumps in discs with and without 
corona, but both with M = 10 s Mq. We see that not only 
the self-gravitating disc decreases because of the corona, but 
also the mass of the clumps increases. But the effect of a 
corona on the instability of the disc is more effective for 
ao = 0.3 comparing to ao = 0.03. It means as the viscos- 
ity coefficient ao increases, the corona makes the disc more 
gravitationally unstable. In fact, in our model, the fraction of 
the dissipated energy into the corona is directly proportional 
to the viscosity coefficient. So, as this parameter increases, 
since more energy dissipated in the corona, we see the effect 
of the corona more effectively. Comparing the cases with a 
corona and without a corona for ao = 0.3, we see that i? ag 
decreases by a factor of approximately 2, but Mf rag increases 
with more or less the same factor because of the existence 
of the corona. 

The above behaviors are also seen in Table 2, but the 
mass of the first clumps are much smaller in the case of 
M — 3 x 1O 6 M0 which corresponds to Sgr A*, comparing to 
the other case. Also, in this case, the self-gravitating radius 
is much larger. However, evidently corona increases the mass 
of the first clumps and reduces the self-gravitating radius. 
Massive stars were found at a distance of order 0.1 — 0.3 pc 
from Sgr A*. For this central back hole, we have 0.1 pc ~ 
3 x 10 Rs- According to Table 2, the self-gravitating radius 
of our models is within the range 0.03 to 2.3 pc, depending 
on the existence of the corona and the input parameters. 
For example, consider a disc without corona with Ze /e = 10 
and Qo = 0.3, this table gives i? sg ~ 0.5 pc. But when we 
consider the same disc but with corona, the self-gravitating 
radius decreases to f? sg ~ 0.25 pc. This new location of the 
first clumps is consistent to typical distance of the observed 
massive stars in the Galactic center. Nayakshin & Cuadra 
(2005) found the self-gravitating radius for their standard 
disc is around 0.03 pc. We note that their standard disc 
model is not similar to our standard model, but since with- 
out a corona, we calculate a larger radial distance for the 
first clumps, if we modify the model Nayakshin & Cuadra 
(2005) to include a corona, then the self-gravitating radius 
decreases. If we consider higher values for the viscosity coef- 
ficient, the discrepancy between the cases with and without 
corona would be more significant. 
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Figure 3. The Toomre parameter for a disc with a corona (solid 
lines) and without a corona (dashed lines) vs. radial location in 
the disc with n = 1.0, Mm = 0.6 and k = 1 for M = 10 s A/ Q (top) 
and M = 3 x 1O 6 M (bottom). 



4 DISCUSSION 

We have studied the effect of corona on the gravitational 
stability of the disc by assuming that the hot corona and 
the disc coexist at the same radius from the central black 
hole. In our disc-corona model almost all the mass exists in 
the disc, but the dissipated energy in the corona depends 
on the input parameters such as mass accretion rate, mass 
of the central black hole and the exponent of the viscosity 
prescription. Our study shows that corona has a destabiliz- 
ing effect on the gravitational stability of the disc. However, 
our disc and corona system is still able to explain the for- 
mation of rings of massive stars observed in the Galactic 
center and in the nucleus of M31 (Davies et al. 2006). But 
star formation in self-gravitating disc with corona may have 
significant differences with discs without corona. Consider- 
ing results of our study, one can try to determine the initial 
mass function for disc-born first clumps at the inner region 
of a self-gravitating disc with corona around a super massive 
black hole, though Nayakshin (2006) has already predicted 
a top-heavy initial mass function for the first starts at the 
distances between 0.03 and 0.3 parsec from our Galactic cen- 
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ter. We think that the estimated supermassive mass of the 
fragments according to Goodman & Tan (2004) will change 
because of the destabilizing effect of corona which transports 
part of the generated energy inside the disc to the corona. 

As we discussed, there are significant observational ev- 
idences in favor of star formation activities in accretion disc 
of AGNs. The mechanism of accretion is an important part 
of any theory for accretion discs. Most of the analytical and 
numerical simulations show that MRI is the most signifi- 
cant process in accretion discs. Almost all models of star 
formation in accretion discs use MRI as the main physical 
process of accretion. Does MRI have other consequences for 
the structure of the disc? As we discussed in the introduction 
section, the idea of a hot corona lying above the AGN accre- 
tion disc has been mainly developed in order to explain sub- 
stantial emission in the X-ray range. There are interesting 
analytical and numerical studies for understanding processes 
of corona formation (e.g., Heyvaerts & Priest 1989; Miller & 
Stone 2000). Interestingly, current models of disc and corona 
systems replace the fraction of accretion power transferred 
from the disc to the corona with a Poynting flux quantity es- 
timated from a mean field buoyant velocity and an equipar- 
tition, mean field magnetic energy density. In other words, 
MRI may lead to corona formation as well, even though all 
the details of the process itself are not yet understood. Hey- 
vaerts & Priest (1989) have shown, loops and arcades can 
form also in AGN accretion discs; these structures, connect- 
ing remote points of the disc itself, can convert disc kinetic 
energy into magnetic energy and subsequently dissipate it 
by emitting the observed spectrum. Numerical simulations 
indicate that turbulent fluctuations in a vertically stratified 
disc are capable of driving the magnetogravitational modes 
of the Parker instability (e.g., Miller & Stone 2000). It seems 
there is a coexistence between MRI and corona, at least in 
AGNs which its X-ray emission can be explained based on 
the existence of corona. Since our model is also based on 
this idea that MRI inside of the disc is also the source for a 
magnetically dominated corona, we think, this model is ap- 
plicable to such AGNs, if the accretion disc itself is massive 
enough. 

One should note that we have not followed the evolution 
of the first clumps after formation. One of the main impor- 
tant factors is accretion onto these new fragments. In fact, 
although the typical mass of the first clumps correspond- 
ing to the Galactic center is low comparing to the cases 
with higher central mass, we think, the accretion onto first 
clumps is able to increase the mass of first clumps signif- 
icantly. This accretion process is believed to be similar to 
the growth of terrestrial planets in a planetesimal disc. In 
order to estimate, the amount of the accreted mass onto 
the first clump, we should calculate the Hill's radius which 
is directly proportional to the mass of the clump itself at 
the time of formation. Since corona increases the mass of 
the first clumps, the corresponding Hill's radius and subse- 
quently the rate of the accretion onto it will increases. It 
means the final mass of the first clumps is even much higher 
because of the corona and subsequent enhanced accretion 
rate comparing to a case without corona. 

Star formation feedback is indeed able to slow down 
disc fragmentation as suggested by several authors, but since 



our goal was only to study the possible effects of the corona 
on the gravitational instability of the disc, we did not in- 
cluded star formation feedback in our model. Also, we have 
not studied the evolution of the clumps after formation. In 
fact, we showed that corona may have significant effects on 
the star formation in AGN disc or even near our Galac- 
tic center which should be considered in any successful the- 
oretical model. Since our model for disc corona system is 
a phenomenological model, more careful description of the 
disc-corona transition is necessary. This transition zone may 
have an important effect on the emergent radiation flux. On 
the other hand, we did not investigate the structure of the 
corona itself, because in our simple model there is only en- 
ergy exchange between corona and the underlying optically 
thick disc. In future works, it would be interesting to study 
gravitational stability of a system of disc and corona, in 
which not only the corona itself can transport the angular 
momentum, but also there are exchange of mass, energy and 
angular momentum between disc and corona. 
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